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I. INTRODUCTION 

It has been well known that spontaneous (continous) 
symmetry breaking states usually support gapless Gold- 
stone modes which give arise various quantum fluctua- 
tions of local order parametersi^*^. In one dimension 
at zero temperature, long wave length quantum fluctua- 
tions can be overwhelming and extremely harmful to the 
establishment of long range order. Depending on details 
of energetics and topology of order parameter manifolds, 
quantum fluctuations might result in complete symme- 
try restoring by driving the states away from the spon- 
taneous symmetry breaking solutions, or lead to alge- 
braically long range order instead of the conventional off- 
diagonal long range order in high-dimensional systems. 

For 0(3) nematic states of spin-one bosons in optical 
lattices (either condensates or Mott states) studied re- 
cently in quite a few articles£i£i£*2*i2iii*i£, indeed in one- 
dimensional lattices the symmetry is completely restored 
and resultant states are rotationally invariant. Unlike 
high-dimensional nematic states, spin excitations in these 
states are fully gapped as a result of symmetry restoring. 

From the point of view of (1 + l)d Euclidean space, 
one can also attribute, at least partially, the disruption 
of long range order to the proliferation of Meron- or 
Skyrmion-like topological instantons. This fact suggests 
that the topology of vacuum manifold (mean field) should 
play an important role in maintaining or destroying long 
range order. 

For instance, when a magnetic field is applied to 0(3) 
nematic states, the easy axis of the nematic order pa- 
rameter tensor is pinned in a plane perpendicular to the 
external field. This lows the symmetry of nematic order 
and the 0(3) symmetry is reduced to 0(2) one, i.e. the 
state becomes a planar nematic state. For nematic con- 
densates this mechanism leads to a vacuum manifold of 
[S 1 x S 1 ]/Z 2 instead of [S 1 x S 2 ]/Z 2 in the absence of 



magnetic fields. 

The above observation then implies that in one- 
dimension in the presence of external magnetic fields long 
wave length fluctuations should not completely disrupt 
nematic ordering in a plane perpendicular to the exter- 
nal fields as pointed out in a previous article^. Instead 
one should expect an algebraically ordered nematic spin 
liquid. In this paper we are going to investigate this pos- 
sibility from a very different starting point and develop a 
theory which fully and quantitatively characterizes field- 
induced algebraically ordered nematic spin liquids in one- 
dimensional optical lattices. 

A closely related subject, magnetically stabilized ne- 
matic order in three-dimensional bipartite lattices, has 
also been studied recently 12 . It was found that the sta- 
bilization of nematic order can be studied via a map- 
ping to the ferromagnetic XX Z model and the critical 
phenomenon is identical to condensation of ideal bosonic 
magnons. Away from critical fields the nematically or- 
dered states are described by weakly interacting dilute 
magnons. 

We are going to employ the previously developed map- 
ping between the low energy Hilbert space close to critical 
fields and the ferromagnetic S — 1/2 XX Z model to fur- 
ther analyze one-dimensional field-induced algebraically 
ordered spin nematic liquids. Though the general ap- 
proach employed below is similar to that developed for 
three-dimensional lattices, the unique feature of dilute 
magnons in the one-dimension limit makes the current 
situation distinct and physics qualitatively different from 
that discussed in the previous article. In fact, it is the 
distinct properties of fermionized dilute magnons which 
lead to anomalous critical exponents as well as fascinat- 
ing algebraically ordered nematic liquids discussed below. 

In Sect. II, we are going to review some basic re- 
sults of the mapping between the field-driven quantum 
phase transitions (between spin singlet Mott states and 
nematic Mott states) and the ferromagnetic S = 1/2 
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XXZ model. We further demonstrate the linear re- 
lation between the pseudo-spin operators in the XXZ 
model and nematic tensor order parameter operators. In 
Sect. Ill, we study the critical behavior of these quan- 
tum phase transitions in the Holstein-Primakov boson 
model, and point out that one-dimensional quantum fluc- 
tuations strongly renormalize the critical exponents. Un- 
like in three-dimensional lattices where the critical point 
is fully characterized by ideal bosons (noninteracting) , 
in one-dimensional lattices the critical point is charac- 
terized by hardcore bosons with an infinitely large inter- 
action constant. And we believe our evaluation of the 
critical exponent is exact. 

In Sec. IV, we further investigate the critical point 
physics using the Jordan- Wigner's fermionic representa- 
tion. We obtain an identical critical exponent as in Sec. 
III. We argue that the critical point is also fully charac- 
terized by the free Fermion model; however the Friedel 
oscillations, one of the basic aspects of Fermions, do not 
appear in the correlation function of nematic order oper- 
ators. 

In Sec. V, we further study the one-dimensional al- 
gebraically ordered nematic spin liquids close to critical 
fields using a strongly interacting boson model. We em- 
ploy the technique of bosonization to evaluate various 
correlation functions in the nematic spin liquids and ar- 
gue that in general these liquids, up to the Friedel oscilla- 
tions are equivalent to the Luttinger liquids of Fermions. 
In Sec. VI, we conclude our studies and discuss some 
consequences of these theoretical findings. 



even number of particles per lattice site, though all con- 
clusions arc qualitatively correct for an odd number of 
particles per site. 

In high-dimensional lattices, when the external field is 
absent and n — J ex /E s is smaller than a critical value n c , 
the ground state is a spin singlet Mott state for an even 
number of particles per site (N). And when r) is larger 
than r/c, the ground state is a nematic Mott state. The 
critical value r\ c has been calculated in various limits and 
turns out to be of order of unity 10 . 

In one-dimensional lattices, in the absence of external 
magnetic fields the ground state is always a spin singlet 
one and the mean field solution with nematic ordering 
is unstable. When J ex is set to be zero, at each lattice 
site the lowest energy state is a spin singlet one (S = 0) 
and the first excited states are fivefold degenerate 5 = 2 
states. The ground state in this case is simply a product 
of spin singlets at each lattice. 

When an external field is applied along the z-axis, one 
of the five-fold degenerate S = 2 states approaches the 
spin singlet and level crossing takes place. Generally 
speaking close to the level crossing points at each site, 
two states \S,S Z > and \S + 2, S z + 2 > become degen- 
erate; the on-site low energy Hilbert space is spanned by 
the two states involved in the level crossing. 

It is therefore convenient to express the effective low 
energy Hamiltonian in this pseudo-spin subspace. Con- 
sider the crossing between states |0, > and |2,2 >. If 
one introduces the on-site pseudo spin Hilbert space Sk 
which consists of two states, 



II. THE EFFECTIVE THEORY CLOSE TO 
CRITICAL POINTS: THE FERROMAGNETIC 
5=1/2 XXZ MODEL 

We introduce ipj. a (tpka), a = x,y,z &s the creation 
(annihilation) operators for a spin-one particle in three 
orthogonal states at site k. The effective Hamiltonian 
for Mott states of spin-one bosons with antiferromagnetic 
interactions in optical lattices in an external field can be 
conveniently expressed asi^^H^, 



| T) = \S = 2, S z = 2), 1 1) = \S = 0, S z = 0), (2) 

then the truncated Hilbert space for the whole lattice St 
is the product of pseudo spin Hilbert space Sk at each 
site k, i.e. St — Si <8> S2 <8> ... ®5t® .... 

Defining | |) and | J.) to be two eigenstates of Pauli- 
matrix <7 Z , o z \ |) = | |), o z \\) = — | j), one obtains the 
effective Hamiltonian which is represented by the ferro- 
magnetic XXZ model in an effective external field along 
the z-directioni^, 



H - E s ^2 Sfe - SkzH z 

fc k 

-Jex^(Qa/3{k)Ql3a(l) + h.C.)\ 
(kl) 



xz 



(1) 



Here E s is the "bare" spin gap studied in Ref.lQ) and 
J ex is the strength of exchange interaction. Sk a — 
—itufi^Xc^kp is the total spin operator for lattice site 
k. Q a p(k) is a bilinear tensor operator, the expectation 
of which yields the nematic order parameter. Finally, 
the sum over (kl) represents a sum over all neighboring 
sites. And again we are interested in the case with an 



{kl) 



0~kaO~la 



f3y^^O~kzO~lz 

(kl) 



h z ^2a kz . (3) 



Eq.Q holds when the truncation is valid, which requires 
that J ex should be much less than E s . Hereeo, (3 and h z 
were calculated in Ref . lfl2T l . 

The Hamiltonian in Eq.© has an 0(2) invariance in 
the xy-plane corresponding to the 0(2) rotational invari- 
ance in the microscopic Hamiltonian in Eq.Q) when the 
external magnetic fields are applied along the z-direction. 
This 0(2) invariance represents the 0(2) nematic sym- 
metry we are going to examine. Furthermore, in general 
the truncation can be applied in the vicinities of all crit- 
ical points where level crossings between \S,S Z = S) and 
\S + 2, S z = S + 2) occur, S + 2 < N. 
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The main issue dealt in this work is the induced spon- 
taneous symmetry breaking in the presence of external 
magnetic helds. Particularly we are interested in the 
breaking of 0(2) rotational symmetry in the xy-plane 
when a field is applied along the z-axis. Unfortunately 
all operators that are linear in ip^, or ip a vanish iden- 
tically in Mott states regardless whether the rotational 
symmetry is broken or not; so linear operators fail to 
distinguish two phases which interest us. 

The rank-two-tensor order parameter introduced on 
the other hand vanishes only in rotationally invariant 
spin singlet states and develops nontrivial structures 
once states break the 0(2) symmetry. For this rea- 
son we introduce the projected nematic order parameter 
Q^p{ot, (3 — x, y, z) which is designed to signify the 0(2) 
symmetry breaking in the xy plane as follows 



L /3'a' 



(4) 



Two constant tensors we would like to project away are 
defined as^i 



n a/3 = ^/|*ea/37 S 7: Ha/3 = ~^ SaS P ~ 3<W)' ( 5 ) 

and s is the unit vector along the direction of spin po- 
larization which in our case is along the z-axis. It is 
evident following the discussions in Ref. |l2T) that the vac- 
uum manifold (mean field) of field-induced spin nematic 
states is S /Z%, identical to that of the standard planar 
nematic order parameter. 

In the pseudo spin subspace, the above operator can 
be expressed in terms of the Pauli matrices a , 



?a/3 



= cr + r a/3 + (7 r^,a,/3 = x,y,z; 



(6) 



and T a f3 is a constant symmetric traceless 3x3 matrix. 
The detailed form of the T-matrix depends on the number 
of particles per site as well as the level crossings involved. 
One can find more discussions on the calculation of the 
T-matrix in Appendix A. 

For instance for 2-particles per lattice site, one finds 



- a/3 



and for four-particles per site, one obtains 




(7) 



1 



-i 



IV = A /- -i -1 
ltM 



(8) 



It is more convenient to further truncate the projected 
nematic tensor operator Q^a m the ^y-plane which in- 
terests us by only keeping the elements with a, (3 = x,y. 
The truncated operator then has the following structure 



APxy 



'-a/3 



(9) 



and T a p is a constant symmetric traceless 2x2 matrix 
obtained by truncating the T-matrix. For 2-particles per 
lattice site, 



a/3 



1 

V3 



1 -i 

-i -1 



(10) 



Finally, the magnetization can also be expressed in 
terms of the pseudo spin operators. For instance, the 
z-component is simply proportional to a z . Following the 
definition of the pseudo spin operators, one finds that 



M fc 



Cfc2 



1. 



(11) 



(Here we have set the Bohr magneton as one). The rela- 
tion between the nematic tensor operator and Pauli ma- 
trices, together with the effective Hamiltonian allow us 
to investigate various nematic spin correlations in terms 
of the XX Z model. 

We also would like to point out that the subspace 
spanned by pseudo spin states is actually a polarization 
(in the xy plane) free subspace. That is 



< S a >— Q, a = x, y. 



(12) 



The absence of polarization (in the xy plane) is a usual 
property of a planar spin nematic state. 



III. CRITICAL STATES AS ID 
HARDCORE-BOSON GASES 

As demonstrated in a previous work, the field-driven 
quantum phase transitions between spin singlet and spin 
nematic Mott states can be studied in terms of the 
condensation of spin-two magnons in spin singlet Mott 
statesi^. To carry out discussions along this line, it 
is rather convenient to introduce the Holstein-Primakov 
representation for the XXZ ferromagnetic spin model. In 
this representation, all spin operators are expressed in 
terms of Holstein-Primakov bosons, 



c'c) c 
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c^V 1 — etc 
2c t c - 1 



(13) 



(c) is the creation (annihilation) operator of bosons sat- 



«tl = 



ct 

isfying the usual bosonic commutation relations 
1; and the raising and lowering operators are defined as 
<j + = (a x + io~ y )/2, o~ = [a x — io~ y )/2. One can further- 
more verify that [o~ a , ap] = i2e a p 1 o~ 1 and a ■ a — 3. The 
Hamiltonian of the XXZ model then transforms intc 12 



4 



H 



xxz 



2 J2 V t 1 ~ C fc C fc] C fc C / V 1 ~ 4 C * 



<fci) 

+ (4(1 + /3)-2fc,) ( 14 ) 



Previous calculations in Ref.jJJj) show that (3 is always 
negative for interacting spin-one bosons under consider- 
ation. 

In the dilute gas limit which interests us, the number 
of Holstein-Primakov bosons per lattice site is much less 
than one i.e. n c — (c^Cfe) <C 1. For this reason, one 
can expand the nonlinear operators of cr ± in terms of n c . 
Indeed, the resultant many-body Hamiltonian up to the 
second order of n c is 



Hxxz 

J ex 



J^[4 - 4cosq + 4/3 - 2h z ]c\c q 



V v b(<i)4i+A 



q 2 -q C <!i C 92J 



V B (q) = -4/3 COS q + 0(p 2 ). 



(15) 



The lattice constant a has been set as one. And at small 
q the dimensionless mass of bosons m is 1/4. 

The Hamiltonian in Ea.lfT5)l. or Ea.lfTl)l in fact de- 
scribes bosons interacting with short range interactions. 
The bare two-body repulsive interaction potential Vs(fcZ) 
when two particles are in two neighboring sites k and 
I — k ± 1 is — 4/3 (in units of eoJ ex , and — 1 < /3 < 0) i.e., 



V B (kl) = -4p6 k±hl . 



(16) 



Therefore, the g-dependent two-body interaction con- 
stant becomes —4/3 as q approaches zero. 

Following Eg. 1151) . when h z < h c — 2/3, all magnon 
excitations are fully gapped with the energy gap being 
2(h c — h z ). The ground state therefore is a vacuum of 
the bosonic field or a spin singlet Mott state; at h z = 
h c = 2/3, excitations with q = become gapless signifying 
a quantum phase transition between spin singlet states 
and nematic states. When h z > h c , magnons appear in 
the ground state and the density is determined by the 
following standard equilibrium condition 



2h, 



2h c = ii{(3,p). 



(17) 



Here p is the chemical potential depending on the density 
p and the parameter (3. Therefore nematic Mott states 



close to critical points are effectively described by dilute 
interacting bosons (spinless). Ea. l(T7|) indicates that close 
to critical points h z — h c approaches zero and bosons are 
in the dilute limit. 

A similar set of equations have been derived and stud- 
ied in Ref.(12). In high-dimensional lattices, one can 
easily apply the results on weakly interacting bosons to 
analyze field-induced nematic Mott states at and close 
to critical points. This has been achieved in a previous 
work and the relation between the development of ne- 
matic order and condensation of magnons has been fully 
investigated^. The chemical potential of a weakly in- 
teracting magnon condensate with interaction given in 
Ea.lfTojl is simply — 8/3p. Following Ea. (|T7J) . this leads to 
the magnon density as well as the magnetization close to 
the critical field. 

However, in one-dimensional lattices with Vt lattice 
sites and a given number of bosons, the single-particle 
zero point motion energy scales as V^ 2 while the inter- 
action energy per particle scales as Vf 1 . As a result, the 
dimensionless coupling constant 7 which is defined as the 
ratio between the interaction energy per particle (— 4p/3) 
and the kinetic energy of a particle confined within a 
length scale p^ 1 (p 2 ) is inversely proportional to p, the 
number of particles per lattice site, 



7 



4/3. 

5 

P 



(18) 



and 7 diverges as p approaches zero in the dilute limit 
close to critical points. One-dimensional dilute bosons 
with interactions arc equivalent to hardcore bosons as the 
density p is taken to be zero and 7 becomes infinite. The 
critical states in one dimension are distinctly different 
from free bosons. 

As pointed out by Girardeau^, the ground state wave 
functions of hardcore bosons (7 — 00) in the continuum 
limit can be expressed as the absolute value of the Slater 
determinant of free fermions. Liniger and Lieb further 
investigated cases when 7 varies from zero to infinity us- 
ing the Bethe Ansatz method^. In addition to confirm- 
ing the solution in the impenetrable hardcore boson limit 
obtained by Girardeau, they also pointed out the impor- 
tant functional relation between the two-body scattering 
phase shift ((71,(72) (see more discussions in Appendix 
B) and the chemical potential of interacting bosons for 
an arbitrary 7. Especially, at large 7 limit the chem- 
ical potential of dilute bosons approaches the value of 
the chemical potential of free fermions with an identical 
density. 

To understand the development of magnetization close 
to critical points, we outline the calculation of the chem- 
ical potential of the Holsctin-Primakov bosons in the di- 
lute limit. We first notice that in the continuum limit 
the two-body interaction in Ea. l|16|l becomes delta-like 
with strength being equal to —4/3. Following Ref^^ (see 
more discussions in Appendix B), we find that the chem- 
ical potential for bosons with the two-body interaction 
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-4(3S(x- 
tions 



■a;') is determined by the following set of equa- 



dqn(q) = p, 
N f m9 2 

-8(3 I dp- 



iip) 



16/3 2 + (p - q) 2 



2im{q) - 1. (19) 



Here n(q) is the distribution function of qausi- 
momentum; the dimensionless mass m is equal to 1/4. 
N is the total number of particles. The first two equa- 
tions above are the definitions of the number density of 
particles p and the total energy E$ of the system. The 
last one in Eq . (|L9I) follows the Bethe Ansatz solution. 
This set of equations have been studied in Refflfll and for 
the convenience of readers we present the derivation in 
Appendix B. The Bethe Ansatz solution for the XXZ 
spin model was also obtained by Yang and YangSS. 

The discussion in Appendix B indicates that the solu- 
tion of Eq. l|19|) should yield the following chemical po- 
tential in the dilute limit 



dE 27r 2 37 3 + 2 7 2 2 
-P ■ 



dN 



3 ( 7 + 2)3 



(20) 



And 7 is defined in Eci. (|18fl . In the presence of external 
fields, in equilibrium the chemical potential has to be 
equal to 2h z — 2h c (> 0) as suggested in Ea. l|17|) . This 
leads to the equilibrium density of strongly interacting 
magnons close to critical points where the dimensionless 
interaction constant 7 diverges as indicated in Eq. l|18|) . 



(21) 



p(h z ) = -\Jh z - h c + ... 
Taking into account Ea. pi(l . one obtains 

M z = 2c{c k , (22) 
and derives the magnetization close to the critical point 



critical exponent found for high-dimensional lattices (See 
Ref|l3). 

Slightly away from critical points, we also obtain the 
corrections to the magnon density as well as the mag- 
netization by expanding Eq.J3U|) in terms of a large but 
finite 7. In the dilute limit, the chemical potential ap- 
proximately is 



2^2 
3 



3P 2 



(24) 



The leading correction to the magnetization M z 
Eq.® is 



SM Z = 



3tt 2 



-(K-h c ) + - 



(25) 



which becomes significant when h z — h c becomes compa- 
rable to (3 2 . 

In Sec. IV, we are going to further examine the crit- 
ical states using the Jordan- Wigner Fermionic represen- 
tation. In Sec. V, we apply the general bosonization 
approach for one-dimensional fluid to the states close to 
the critical points. Finally we would like to mention 
that both bosonic and Majorana fermionic representa- 
tion have been employed to study antiferromagnetic spin 
chains in an external fieldi 3 ^^*^*^. Bethe Ansatz solu- 
tions of the XXZ model have also been used to calculate 
critical exponents in various models^. 



IV. CRITICAL STATES AS FREE-FERMION 
GASES 

Since the ground state wave functions of interact- 
ing bosons in the dilute limit, or equivalently hardcore 
bosons can be related to the Slater determinant of free 
Fermions, not surprisingly many aspects of the fixed 
point solutions reminisce free Fermions. It is certainly 
appealing to work with the fermionic representation (of 
Jordan- Wigner type) and find out to what extend the 
critical points can be characterized by the free Fermion 
model. In this section we are going to implement this 
idea. 

In the Jordan- Wigner representation, spin-operators 
are expressed in terms of string operators defined by 
Fermions, 



M z = -y/h z -h c . 

7T 



(23) 



We should emphasize here that because of the large 
dimensionless coupling constant of dilute bosons in one- 
dimensional lattices, the chemical potential is a quadratic 
function of the density of bosons instead of a linear func- 
tion as in the weakly interacting limit. As a result, 
in one-dimensional lattices the critical exponent in the 
quantum phase transitions between the nematic and spin 
singlet Mott states is renormalized to be one-half of the 



4=exp(m^///0/2 

Kk 

a k = /feexp(-i7r^///0 



Kk 



Okz 



Vlfk 



1. 



(26) 



The string operator in the Wigner- Jordan form is ori- 
ented along the path starting at —00 and ending at k—1. 

Substituting the expressions for spin operators, one 
rewrites the ferromagnetic XXZ model in the following 
form (see more discussions in Appendix C) 



6 



H 



xxz 



eJ P , 



-2 £ (flh + /// fc ) 



<fei> 



+ [4(l + f3)~2h z ]J2flfk 
k 

-4(1 + /3) £ flfkfifi- 



(27) 



<ki> 



Eg. 1)27(1 shows that the ferromagnetic XXZ model 
(0 > /3 > — 1) is equivalent to attractively interacting 
Fermions. The bare two-particle attractive interaction 
Vp(kl) is short ranged; the interaction is nonzero only 
when two particles occupy two neighboring sites and has 
a negative value of —4(1 + (3), i.e. 



V F (kl) = -4(1+0)5. 



k±l,l- 



(28) 



Obviously the bare potential vanishes only when /3 = — 1. 

However, because of the Pauli exclusion principle, it is 
well known that scattering between spinless Fermions is 
absent if the two-particle potential is truly S- like. For a 
similar reason in a dilute Fermi gas the scattering is very 
weak if the two-body potential is short ranged. There- 
fore, in one-dimension for a given number of Fermions 
while the one-particle kinetic energy for Fermions is pro- 
portional to Vff 2 as in the case of bosons, the interaction 
energy per Fermion is proportional to Vrj7 3 instead of Vf 1 
as a result of the Pauli exclusion principle. It indicates 
that the dilute limit of Fermions is truly weakly interact- 
ing and the fixed point of free Fermions is stable. Indeed, 
in one-dimensional lattices, the dimcnsionlcss coupling 
constant scales as pa in the dilute limit (a is the range of 
interaction which in our case is the lattice constant; it is 
set as one in this article). 

Following the discussions in Appendix C, one obtains 
the following effective Hamiltonian 



H 



xxz 



= ^[4/3 - 2/i z + 4 -4cos <?]/]/, 



V F {q) = 4(l + /3)(cos?- 1) <0 



(29) 



And again at the critical point h c = 2/3, the excitations 
with q = become unstable and the quantum phase tran- 
sition occurs. When h z < h c , Fermionic excitations are 
gapped with an energy 2(h c — h z ) and the ground state 
is a vacuum of Fermions representing the spin singlet 
Mott state. Slightly above the critical field, Fermions 
with small momenta have negative energies and we have 
a dilute interacting Fermi gas in the ground state repre- 
senting nematic Mott states. 

The second sum in Ea. (|29|) indicates that the forward 
scattering amplitude between right-movers or left-movers 



in one-dimension, i.e. (rqF, rq F ) — > (rqp^rqp) (r = ±) 
is proportional to q 2 ~ (q denotes the transferred mo- 
mentum during the scattering and qp is the Fermi mo- 
mentum) and therefore vanishes in the long wave length 
limit. However the back scattering between right-movers 
and left-movers, (rqp,—rqp) — > (—rqp^rqp) involves a 
momentum transfer q = 2qp\ the scattering amplitude 
is proportional to qp and remains finite in the long wave 
length limit. As expected before, the dimensionaless cou- 
pling constant which is the ratio between the interaction 
energy e(l + (3)p and the kinetic energy e should be pro- 
portional to 



lF = -(l + (3)p. 



(30) 



Notice that the bare two-particle attractive interac- 
tions —4(1 + (3) are nonzero for any system with (3 =/= — 1 
and can even be comparable to the Fermi energy. In 
principle, Fermions under consideration can be in the 
strongly attracting limit. However, the effective coupling 
constant which itself depends on the density of Fermions 
is always small in the extremely dilute limit. And it be- 
comes precisely zero when critical points are approached 
and one anticipates that the Fermion density vanishes. 
Close to the critical points, we do have a collection of 
weakly interacting Fermions, a problem which can be eas- 
ily handled. Of course when /3+1 = 0, the bare two-body 
interactions become zero identically. 

To summarize the result of this mapping, we would like 
to emphasize that all states along the triangular bound- 
aries in the phase diagram in Fig.l a) correspond to the 
free Fermions. The 0(3) symmetric tricritical point on 
the other hand is a free boson critical point. All states 
within the triangular boundaries are Luttinger liquids 
with varying densities and Luttinger liquid parameters. 
Finally states in the SSMI phase outside the triangular 
boundaries in Fig. lc) are represented as the vacuums of 
Jordan- Wigner Fermions and the SPMI corresponds to a 
filling factor N — 1 Fermionic Mott state (see Fig.l for 
more discussions). 

Therefore, the critical point can also be studied conve- 
niently in the Fermionic representation. For instance, the 
chemical potential in this case is naturally proportional 
to the square of the Fermion density which is consistent 
with the results obtained in the bosonic representation. 
Taking into account that m = 1/4, one obtains 



Pf = 2ir 2 p 2 



(31) 



Along with the equilibrium condition for Fermi gases 
(which is identical to the condition for bosons in Eq.l|17|)l 



Hf{/3,p) = 2(h z - h c ), 



(32) 



one obtains the Fermi density close to the critical point. 

Ea. i|ll[l and Ea. H2tj|) further indicate that the magne- 
tization is proportional to the Fermion number operator 
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M kz = 2flf k . 



(33) 



This again leads to Eq.(P3J. 

Moreover, one can express nematic tensor operator in 
terms of the Fermion operators 



Qlf = T aP exp(*7r flfOfl + h- 



(34) 



Kk 



In deriving it one has to take into account Eq.^and 
Ea. (|26JI . the mapping between the pseudo-spin operator 
and Fermion operators. Note the nematic operator has 
nonlocal dependence on the Fermion operators: it is a 
string operator of Fermions! This aspect of the tensor 
operator plays a very important role in the nematic cor- 
relations in ID. 

So we can use the free Fermion properties to calculate 
the nematic tensor correlators in ID. Indeed, we find that 
the correlator appears to be clearly related to the free 
Fermion correlator, 

<ori j/ wQ^(o>= 

< fk exp(-i7T ^ flfn) exp(l7T ^ fmfm)fi > +CC. 



rn </ 



(35) 



Evaluating this correlator is tricky because a Fermion 
operator itself inexplicitly contains string operators and 
is correlated with the string operator. The consideration 
towards the end of Sec. V and in Appendix C implies that 
the contribution from the string operator in the Jordan- 
Wigner representation yields precisely a sign oscillation 
in the Fermion wave function. So the appearance of the 
Jordan- Wigner phase factor in the nematic correlator en- 
sures that such oscillations in the free Fermion correlator 
get cancelled out. Finally, standard calculations yield 



i 



|iA)o 



(36) 



where the Friedel oscillations that appear in the free 
Fermion correlator are indeed absent in the nematic cor- 
relations (We have also neglected the higher harmonics) . 
?/o = 2 following the Luttinger liquid theory (see Ap- 
pendix D). In the next section, we are going to study the 
nematic states close to critical points using the bosoniza- 
tion technique. 



V. ALGEBRAICALLY ORDERED SPIN 
NEMATIC LIQUIDS CLOSE TO CRITICAL 
POINTS 

Following the discussion at the end of section IV and 
Appendix C, as far as the nematic correlations are con- 
cerned, we notice that the Fermion representation of 



m 









\. MI 






LL 




v 








/ VAC 







p 



m 




FIG. 1: (Color online) a) The phase diagram for the Jordan- 
Wigner Fermions in the h z -f3 plane (/3 £ [— 1, oo]). Inside the 
triangle are the Luttinger liquids (LL). Outside the triangle 
and h > 0[< 0] is the Mott insulator (MI) [vacuum (VAC)] 
of spinless Fermions. Along the boundary of the triangle is 
the free Fermi gas. b) The phase diagram for the Holstein- 
Primakov bosons (ID) in the h z -/3 plane. Inside the triangle 
is the superfluid (SF). Outside the triangle and h > 0[< 0] 
is the Mott insulator (MI) [vacuum (VAC)] of bosons. Along 
the boundary of the triangle is the hardcore boson gas. c) 
The phase diagram for spin-one bosons in ID optical lattices 
indicated by a) and b). Inside the triangle is the algebraically 
ordered nematic spin liquid (AONSL). Outside the triangle 
and h > 0[< 0] is the spin polarized Mott insulator (SPMI) 
[spin singlet Mott insulator (SSMI)]. The value of for spin- 
one bosons was calculated microscopically in RefH^ and turns 
out to be always negative. 
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pseudo spin operators is redundant: the phase string ap- 
pears twice in the calculation. Taking this point of view, 
we find it is more convenient and natural to employ the 
bosonic representation (see section III) for spin opera- 
tors. And we choose to use the bosonization technique 
to evaluate correlations at (h z — h c ) and close to critical 
points (h z > h c ). 

The bosonization scheme is a generic representation to 
address interacting bosons as well as Fermions in one- 
dimension systems. It is particularly powerful when de- 
scribing algebraically ordered quantum spin nematic liq- 
uids close to the critical points. And it also serves very 
well as a framework to study the equivalence between the 
interacting boson model in section III and the Fermion 
(free and interacting) models in section IV; in fact the 
dynamics of algebraically ordered nematic spin liquids is 
identical to Luttinger liquid ones. We are going to work 
in the continuum limit. 

First we rewrite the Hamiltonian in Ea. (|15fl in terms 
of c'(x), c(x) in the continuum limit ( we have chosen the 
lattice constant a to be one and x is dimensionless), 



xz 



dxc 1 (x)(2h c - 2h z 



2m 



dl)c{x) 



4/3 / dxdx'v B (x - xV (x)c(x)J (x')c(x'). (37) 



Here —AI3vb{x—x') represents the short ranged two-body 
potential in the lattice model and in the continuum limit 
it is taken to be a S- function, i.e., vb(x — x') = 8(x — x'). 
The dimensionless mass of bosons is set to be m = 1/4. 

Following Refl2ll we introduce the following standard 
conjugating bosonic fields for the creation operator of 
bosons c*>(x) (See details in appendix D), 



t(a:) = J^J~^^ e i^) £ e 2 ™ e ^(38) 



m— — oo 



[d x Q(x),$(x')] =i-K%5{x-x'). 



(39) 



And p(h z ) has been given in Eq.J2U in section III. 

Substituting this Ansatz into Ea. (|14fl and taking the 
long wave length limit, we find the following Hamiltonian 
in terms of <d(x) and 4>(x) (h c > h z ), 



^S2£L = {2h c -2h z +p)- 

eJex 7T 



dxd x Q(x) 



dx 



1 F)n h 2 
— ^-{d x Q{x)f + — p{hz ){dMx)) 2 



(.40) 



For an unrenormalized theory, the bare value of 
p((3,h z ) is taken to be ~8f3p(h z ) and dp/dp = —8/3. 
Since the theory is renormalizable at the long wave length 
limit, we expect that the low energy effective Hamilto- 
nian should be of the same form as in the bare case except 
that quantities such as p can be strongly renormalized by 



short wave length interactions. In fact this form of the 
effective Hamiltonian is suggested by the hydrodynam- 
ics of bosonic gase a 21 ' 22 , though the 7-dependence of the 
renormalized p and dp/ dp can be obtained only when 
certain microscopies are taken into account. 

Therefore, the effective Hamiltonian only depends on 
the density and the renormalized functional form of the 
chemical potential p{p). And both the chemical potential 
and the compressibility dp/ dp are functions of the dimen- 
sionless interaction constant j((3,h z ) defined in Ea. (|18H . 
Following Ea. (|21[) . one finds that 



47T/3 



Vhz - h c 



(41) 



as a function of h z . 

In equilibrium, the first term in Ea . (j4(J|l which is lin- 
ear in 0(a:) has to vanish because of the stability of 
the system. This sets the known relation between the 
chemical potential and external field, 2(h z — h c ) = p(p) 
which has been discussed in the previous sections. For 
the weakly interacting systems, in the leading order of 7, 
the chemical potential takes the bare value of — 8/3/?; for 
the strongly interacting case as in the ID dilute limit, p 
is renormalized and discussions in section III show that 
p is proportional to 2n 2 p 2 in the leading order of 1/7. 

Taking this into account, we substitute into Ea. H4U|) 
the renormalized chemical potential p. We recast it in a 
more familiar form 



H 



XX z 



— I dx[v N (8 x Q(x)) 2 +vj{d^{x)f 



Here 



v N (/3, h z ) = 



1 dp(p) 
nh dp(h z ) 



Vj = irh 



p{h z ) 



(42) 



(43) 



again m = 1/4. An alternative way of calculating can 
be found in an early workSl. vj here is taken to be inde- 
pendent of interactions as a result of Galilean invariancc 
(approximate) in the continuum limit; in general when 
periodical or random potentials are present, vj could be 
renormalized away from this value. 

To calculate vjv and vj, one has to take into account 
the Bethe Ansatz result for the chemical potential in 
Ea. (|20() . One obtains the following results for vn, 



v N = 4irhp(h z )- 



(-4/3) 4 



-ip + 2p{h z )Y 



(44) 



Eq. (|21|) shows that both vn and vj vanish as h z ap- 
proaches the critical field h c . The discussions in Ap- 
pendix D show that the nematic operator is algebraically 
correlated, and close to the critical point it takes the form 



(q5 9 (*)q5S'(i)} 



TgpToSjyj^ - P)P 

{ P \k-i\yh ■ 



(45) 
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in the long- wave length limit, where the correlation ex- 
ponent is determined by the value of r\ 



7 ? (/3,/ l2 ) = 2 A /^ = 2 
v N 



(46) 



In Fig.l one can find more discussions on various phases 
in the bosonic representation and implications on spin- 
one bosons in one-dimensional lattices. 

The sound velocity of the Holstcin-Primakov bosons 
corresponds to the spin-wave velocity in nematic states. 
Close to the critical points when p(h z ) vanishes, vn ~ vj 
and the spin wave velocity v s = y/vjyvj is 



v s (h z ) ~ <inhp(h z ) = 4ft V h z - h c + 



(47) 



following Eq.J2IJ. 

Ea. p^ and Eo. lfl^ also define the Luttinger liquid 
dynamics of interacting Fermions 23 - 24 - 25,26 . In the stan- 
dard formulation of the Luttinger liquid theory2&, the 
left-moving and right moving Fermion operators are ex- 
pressed in terms of bosonic fields of Q(x) and $(2;) 



fl R = exp(±iq F x) exp(iO(x) ± i$(x)). (48) 



We have also calculated the nematic spin correlations 
using the Bosonization scheme and obtained Zeeman- 
field-dependent algebraically ordered nematic spin liq- 
uids. At critical points, the nematic spin correlation 
functions are identical to two-point correlations of one- 
dimensional free Fermions with the Friedel oscillations 
completely suppressed. Close to critical points, the al- 
gebraically nematically ordered states are equivalent to 
Luttinger liquids of weakly attractive spinless Fermions 
or fluid of strongly repulsive bosons. 

F.Z. thanks CASTU, Tsinghua University, Beijing for 
its hospitality during his stay in the summer of year 2004; 
he would like to acknowledge useful discussions with I. 
Affleck. This project is in part supported by NSF, P. 
R. China(Grant No. 10247002 and 10404015) (H.Z.) and 
in part supported by a grant from UBC and NSERC, 
Canada (F.Z.). 

Note added: The transition described here might be 
considered as an analogue of commensurate- incommen- 
surate transitions and is of quantum Pokrovsky- Talapov 
typeSI. 



APPENDIX A: NEMATIC TENSOR ORDER 
PARAMETER IN THE PSEUDO SPIN 
SUBSPACE 



The dynamics of &(x) and <!>(x) is given by 



H 



L.L. 



^/^ (W) ) 2 + ^n 2 (,)], 



[n(x),$(x')] =ih8(x-x'). 



(49) 
(50) 



Here H(x) = ^dO(x). One can easily verify that Ea. lfl^ 
and Eq.JSHJ) are identical to Eo. lf4*9")) and Eq.JSOJ if one 

identifies that 



n 



2K 




(51) 



This intimate relation between ID bosonic superfluid and 
Luttinger liquids was appreciated before^ 1 .. 



Close to critical points, the Hilbert space is truncated 
into a two-level subspace at each site. Therefore we can 
define a pseudo-spin 1/2 at each site, where | \) denotes 
the state \S + 2, S z = S + 2) and | J.) denotes the state 
I S, S z = S). In this truncated subspace the nematic op- 
erator 



1 



Qa/3 = Ip^l/lp ~ ^S a/ 3i>\lp 7 , 



(Al) 



can be expressed in terms of pseudo-spin operators. 

Consider the case of 2-particle per site. After some 
straightforward calculations, the matrix elements of the 
nematic operator can be obtained as 



<T \Q aP \ T> = 



k i 

-i i 

oo-l 



(A2) 




(T \Q a p\ i) = -= \ -% -1 = (I \Q a p\ T>*, (A3) 



VI. CONCLUSIONS AND DISCUSSIONS 

In this article, we have studied the critical phenomenon 
in the field-driven quantum phase transitions between the and 
nematic Mott state and spin singlet Mott state in one- 
dimensional lattices. The critical states in one-dimension H \Q a a\ I) =0. (A4) 
are in the universality class of impenetrable bosons or 

free Fermions. We have found that because of hardcore Therefore, for an arbitrary state |f2) in this subspace 
nature of magnons in one-dimension the critical exponent 

of magnetization is renormalized to be one-half of the ir>\ — ^ ( • < ^Al t\ ■ ^ ( mC ^\ \ \\ ( 
value in three-dimensional lattices. ' ' — C0S 2 6XP I * 2 / Sm 2 6XP \ 2 J 
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the expectation value of the nematic operator is 

| % o 

(n|Q^,|n) = cos 2 ^ I -% | o 



2 1 3 
- 1 -| 



sin 9 exp (ic, 



1 -i 
_ . -i -1 
2 v 3 V 

1 i 



sin0exp(-i(?!>)— ^= I i -1 | (A6) 
2 v3 \ o 

The first term in the right-hand side of above equation 
is contributed from the fully polarized state and should 
be projected away. The remained two terms give the ex- 
pectation value of the projected nematic operator, which 
is 



(Sl\Q? \n) = sinO exp (i 



sin 9 exp (- 




(A7) 



To express the projected nematic operator in terms 
of pseudo-spin operators, one should look for a suitable 
combination of the tr-matrices which yields the same ex- 
pectation value as Ea. (|A3|) . or Ea. (|A7ll . It is easily to 
verify that Q^p should be represented by the following 
operator in the pseudo-spin space, 

^(ax+Vy)T al 3 + ^(cT x ~i<J y )r* al3 , (A8) 

where 




(A9) 



Dividing two successive equations and equating expo- 
nents, one obtains 

N 

{q l+1 - Qi )L = ^2(e jti - 9 3 . 1+1 ) + 2nl u (B3) 
i=i 

where a set of integers Ii(i = 1, • • •, N) are the quantum 
numbers for this system, and they are all equal to one for 
the ground state. Taking the continuum limit of Eq.{B3j, 
one has 



dp 



n(p) 



lep 2 + (p- q y 



= 2irn(q) - 1, (B4) 



where A is the cut-off momentum which should be de- 
termined self-consistently. Here n(q) is the distribution 
function of quasi-momentum, which satisfies the con- 
straint 



dqn(q) = p. 



(B5) 



And the ground state energy is given by 

A 

N f a 2 
Eo = ~ n{q)^dq. (B6) 
p J 2m 



-A 



Now rescale all qualities in terms of the cut-off mo- 
mentum A, i.e. —4/3 = AA, p — Ay and q = Ax, and 
denote n(q) as g(x). The above three equations take the 
following form 



1 + 2A 



g(y) 



X 2 + (y-x)' 



; dy = 27rg(y), (B7) 



APPENDIX B: THE CHEMICAL POTENTIAL OF 
INTERACTING BOSONS IN ONE-DIMENSION 

The one-dimensional interacting bosons problem can 
be exactly solved by the Bethe Ansatz method. In this 
appendix we will briefly review the Betha Ansatz method 
used in Refill. Let {qi}{i — 1, • • -,N) be the quasi- 
momenta of the N bosons, the Bethe Ansatz solution 
indicates that all the allowed values of qi 's should satisfy 
the following set of equations: 



N 



(-l) JV - 1 e - < * J 



exp i 



Qj) 



(Bl) 



where 9ij is the phase shift. In our case, one finds that 



0fe - Qj) 



-2 arctan 



g* ~ 1j 
-4(3 



(B2) 



and 



7 / g(x)dx = A, 
-1 



#o = 7^-Np 2 ^ I g(x)x 2 dx. 
2m A J 



(B8) 



(B9) 



In the limit of large 7, A is large compared to \x — y\, 
one can then neglect the x dependence of the integrand 
in Ea. l|B7(l : and g(x) turns out to be a constant. The 
integral equation can be solved approximately and the 
result is 



A 

2ttA - 4 ' 



(BIO) 
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Substituting it back to Ea. (|B8|) one obtains the relation 
between A and 7 as 



A = 



7 + 2 



(Bll) 



And finally, using Ea. (|B9|) one can derive the ground 
state energy 

e ° = w ' 2 £(tti) 2 - (B12) 

Consequently, a relation between chemical potentials and 
particle densities is established in this limit, 

dE Q ^ 3 7 3 + 2 7 2 2 
M= -dN = 6^7tT^F P - (B13) 



APPENDIX C: THE JORDAN- WIGNER 
FERMION REPRESENTATION 



1. The effective Hamiltonian for Jordan- Wigner 
Fermions in the XXZ model 



In this appendix we will briefly review the Jordan- 
Wigner Fermionic representation of the XXZ model, the 
Hamiltonian of which is 



2. Calculation of Eq.ll3^1) 

In the continuum limit, can be written in the fol- 
lowing form 



a 
a~ 



exp(~iqFX — iQ(x))p (x); 

f(x) exp(iqpx + i<d(x)). (C6) 



Taking into account the expression for the left- and right- 
movers in Eci. (|48[) . one obtains 

cr + ~ exp(i&(x)) + exp(—i2qpx — i&(x)); 

o~ — exp(-i$(x)) + exp(i2q F x + i$(x)). (C7) 

Therefore the leading contribution to the lowest har- 
monics in the correlation function is 



< cr+er < exp(i$(x) - i$(0)) > 



(C8) 



Hxxz 



-2 E K + a k °i 



<kl> 



-(/? + 1) a kz <7 lz - h z ^ a z 



(CI) 



Following the Jordan- Wigner transformation Ea. (|26|l . 
the spin flipping term becomes 



Calculation of this correlator is the same as those in Ap- 
pendix D. For free Fermions, one obtains the result in 
Eq.(|3BI). 



APPENDIX D: BOSONIZATION AND NEMATIC 
SPIN CORRELATIONS 

In this appendix we briefly review the bosoniza- 
tion procedure for one-dimensional interacting bosonic 
system^. First we express the bosonic operator as 



c'(x) = \J p T {x)e 



(Dl) 



(C2) 



(C3) 



, _ iT E flfl t -*t E ///( 

<°k+i= e i=1 /fc/fc+ie !=1 
= fl(l-2flf k )f k+1 =flf k+1 , 

and similarly, 

a k+i a k = fk+ifk- 

These yield a kinetic term of these Fermions. Neglecting 
a constant term, the second term in Ea. (|Clfl turns out 
to be 

E ak * at * = 4 E fiMlfi - 4 E ( C4 ) 

Thus one arrives at the Jordan- Wigner Fermion repre- 
sentation of the XXZ model as shown in Eq. (|27|) . 

Finally, due to the Pauli exclusive principle one can 
treat the interacting terms as follows 



The mean ground state density is denoted by po, we in- 
troduce a field II(x) to account for the long-wave length 
(> Pq ) density fluctuations. Furthermore in order to 
take the discreteness of particles into account, we intro- 
duce an auxiliary field 0(x) which is related to H(x) via 



-d x Q(x) = Po + IL(x) 



(D2) 



where O increases monotonically by it each time when 
x passes the location of a particle. Therefore the total 
density can be expressed as 



+00 

Pt(x) = -d x e(x) E 6(Q(x)-n*), (D3) 

n— — 00 



where the discrete character is incorporated. Since the 
square root of a 5 function is also proportional to a 6 
function, it yields that 



^2fkfk+lfkfk+l — E yfkfk+lfkfk+l - flflfkf, 
k k 

= E fl-jL+JMcosq-l), (C5) 



+ OO 



(x) ~ Vpo + n(x) E 5{®{x) - nv)e l ^ x 



9,91,92 

and this results in Ea. (|29[l . 



n— — 00 
+00 



Vpo + n(x) 



e 2mie(x) e i<S>(x) 



(D4) 
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where the second equality follows from the Poison's sum- 
mation identity 



+00 



~r^j +OC poo 

E /(»)= E / dzf(z)e 2m ™. (D5) 



m— — 00 



It can be verified that the field and <f> satisfy the com- 
mutative relation 



[d x O(x),$(x')] = mM{x-x'). 



(D6) 



With these expressions, the effective Luttinger Hamil- 
tonian can be read as 

HcS= ^J dx [ v AW x )f + Mdx&{x)) 2 }. (D7) 

Here we have introduced two phenomenological param- 
eters vj and vn- vj is the phase stiffness which equals 
to tmplm in a uniform system, and vn is the density 
stiffness which is defined as 



1 (dp 



nh \dp 



(D8) 



These parameters should be determined by microscopic 
calculations. For boson gases with delta-like interactions, 
these parameters can be obtained from the Bethe Ansatz 
solutions. 

Furthermore we can express the correlation of nematic 
operator in terms of correlators of Holstcin-Primakov 
bosons, 



{c\k)^l - rf(k)c(k)y/l - ct (l)c(l)c\l)) + h.c. 
~ T aP Tl, al (1-p) (ct (fc)c(O) + h.c. (D9) 



The correlation function for the Holstcin-Primakov 
bosons can be computed using the effective Luttinger 
Hamiltonian, 



(J(x)c(0))=p 
+00 , 1 

+ E bm -n 

_ n VPFI 



1 l/(2if) 



p\x\ 

2m 2 K 



{bo 



cos(2irmpx) } . (D10) 



Here 



V V N 



(Dll) 



is a universal parameter only depending on long-wave 
length physics, while the 6 m 's depend on the microscopic 
details. 
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